Conductivity of suspended and non-suspended graphene at finite gate voltage 
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We compute the DC and the optical conductivity of graphene for finite values of the chemical potential by 
taking into account the effect of disorder, due to mid-gap states (unitary scatterers) and charged impurities, and 
the effect of both optical and acoustic phonons. The disorder due to mid-gap states is treated in the coherent 
potential approximation (CPA, a self-consistent approach based on the Dyson equation), whereas that due to 
charged impurities is also treated via the Dyson equation, with the self-energy computed using second order 
perturbation theory. The effect of the phonons is also included via the Dyson equation, with the self energy 
computed using first order perturbation theory. The self-energy due to phonons is computed both using the bare 
electronic Green's function and the full electronic Green's function, although we show that the effect of disor- 
der on the phonon-propagator is negligible. Our results are in qualitative agreement with recent experiments. 
Quantitative agreement could be obtained if one assumes water molelcules under the graphene substrate. We 
also comment on the electron-hole asymmetry observed in the DC conductivity of suspended graphene. 

PACS numbers: 73.20.Hb,81.05.Uw,73.20.-r, 73.23.-b 



I. INTRODUCTION 

The isolation of a single carbon layer via micro mechanical 
cleavage has triggered immense research activity. 1 Apart from 
the anomalous quantum Hall effect due to chiral Dirac-like 
quasi-particlesri 3 - the finite "universal" DC conductivity at the 
neutrality point attracted major attention. 4 For recent reviews 
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The electronic properties of graphene are characterized 
by two nonequivalent Fermi-surfaces around the K and K'- 
points, respectively, which shrink to two points at the neu- 
trality point /i = (/i is chemical potential). The spec- 
trum around these two points is given by an (almost) isotropic 
energy dispersion E(k) = ztvphk with the Fermi velocity 
vf — 10 6 m/s. 11 Graphene can thus be described by an effec- 
tive (2+l)-dimensionalrelativistic field theory with the veloc- 
ity of light c replaced by the Fermi velocity vf^ 

Relativistic field theories in (2+1) dimensions were inves- 
tigated long before the actual discovery of graphen o 13 i 14 and 
also the two values of the universal conductivities of a clean 
system at the neutrality point depending on whether one in- 
cludes a broadening r — > or not were reported then i 15 ' 16 ' 17 In 
the first case, one obtains CpZ^ = §e 2 //ir^ the second case 
yields <7p=Q — ^e 2 /hJ£ Interestingly, the first value is also 
obtained without the limit r — » within the self-consistent 
coherent potential approximation (CPA). 20 We also note that 
the constant conductivity holds for zero temperature, only; for 
finite temperature the DC conductivity is zero^ 

If leads are attached to the graphene sample, an exter- 
nal broadening is introduced and the conductivity is given 
by o-p^ j 22,23,24 which has been experimentally verified for 
samples with large aspect ratio. 25 This is in contrast to mea- 
surements of the optical conductivity, where leads are ab- 
sent and a finite energy scale given by the frequency uj of 
the incoming beam renders the intrinsic disorder negligible, 
T/huj w 0. One thus expects the universal conductivity to be 
given by a^Z® — %e 2 /h, which was measured in various ex- 
periments in graphene on a S1O2, 26 SiC-substrate 27 and free 



hanging. Also in graphene bilayer and multilayers, 2 ^ 8 - as 
well as in graphite 29 the conductivity per plane is of the order 
of er = § e 2 /h. 

The above results were obtained from the Kubo or Laun- 
dauer formula and assumed coherent transport. Also dif- 
fusive models based on the semi-classical Boltzmann ap- 
proach yield a finite DC conductivity at the neutrality point. 
Nevertheless, the finite conductivity was found to be non- 
universa l 30 ' 31 ' 32 ' 33 ' 34 in contraditions to the findings of early 
experiments, which suggested er m i„ « 4e 2 jhX We should 
however stress that one can still assume a certain degree of 
universality, since the experimental values for the conductiv- 
ity are all of the order of 4e 2 /h. It was argued that electron- 
hole puddles 35 or potential fluctuations in the substrate 3 ^ can 
account for a finite conductivity at the Dirac point. An al- 
ternative explanation of this quasi-universal behavior seen in 
experiments is that there is only a logarithmic dependence on 
the impurity concentration due to mid-gap states and therefore 
only in cleaner samples deviations from the universal value 
are seen^ 

On the other hand, the optical conductivity is given by the 
universal conductivity gq for frequencies larger than twice the 
chemical potential fi. It is remarkable that this universal value 
also holds in the optical frequency ranger^ 3 - 6 - a result with 
important consequences in applications ! 37 ' 38 ' 39 Only for fre- 
quencies fuuj < 2/i, the sample-dependent scattering behavior 
of the electrons becomes important and recent experiments in 
show an decay of the universal conductivity with unusual large 
broadening around 2/i which can not be explained by thermal 
effects. 26 Moreover, the spectral weight for ksT <C hut <C 2/i 
does not reach zero as would be expected due to Pauli block- 
ing, but assumes an almost constant plateau of a as uq /3 for 
larger gate voltage. 

The first calculations of the optical conductivity of 
graphene, using the Dirac Hamiltonian were done in Ref. 
H20fl . This study was subsequently revisited a number of 
times ) 41 ' 42 ' 43 and summarized in Ref. lfl7ll . In these calcula- 
tions the effect of disorder was treated in a phenomenological 
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manner, by broadening the delta functions into Lorentzians 
characterized by constant width T. As shown in Ref. 112011 
however, the momentum states are non-uniformly broadened, 
with the states close to the Dirac point being much more af- 
fected by the impurities than those far away from that point. 
In the clean limit, the exact calculation of the optical prop- 
erties of graphene was considered in Ref. a calculation 
recently generalized to the calculation of the optical properties 
of graphene antidot lattices^ 5 . 

In this paper, we generalize the results of Ref. 112011 by 
considering a finite chemical potential, including the effect 
of charge impurities, and the scattering by phonons. We dis- 
cuss two main corrections to the clean system and calculate 
the optical conductivity. First, we include the coupling of the 
Dirac fermions to in-plane phonons, acoustical as well as op- 
tical ones. Out-of-plane phonons only have a negligible ef- 
fect on the electronic properties of graphene. 46 Secondly, we 
include various types of disorder which give rise to mid-gap 
states as well as Coulomb scatterers. 

In Sec. [n] we define the phonon Hamiltonian, deduce the 
electron-phonon interaction and calculate the electronic self- 
energy. In Sec. MI Al we discuss the Green's function which 
is modified due to impurities and phonons. We then present 
our results for DC and optical conductivity and compare it 
to the experiment of Ref. J26ll . We close with remarks and 
conclusions. 



II. ELECTRONS AND PHONONS 

A. Tight-binding Hamiltonian and current operator 

The Hamiltonian, in tight binding form, for electrons in 
graphene is written as 

H = -J2^2t{R,R + 5)[aUR)b<,(R + 5) + H.cll) 

R,a S 

where the operator aj.(i?) creates an electron in the carbon 
atoms of sub-lattice A, and b^R) does the same in sub-lattice 
B. The hopping parameter, t(R, R + 6), depends on the rel- 
ative position of the carbon atoms both due to the presence of 
a vector potential A(t) and due to the vibration of the carbon 



atoms. The vectors 8 have the form 

5i = | (l,V3) <J 2 = |(l,-V3) * 8 = -o(l,0), 

(2) 

where a is the carbon-carbon distance. In order to obtain the 
current operator we write the hopping parameter as 

t -> t e i{e,h)A{tyS . (3) 

Expanding the exponential up to second order in the vector 
potential A(t),and assuming that the electric field is oriented 
along the x direction, the current operator is obtained from 

i (4) 
leading to j x = j£ + A x (t)j^r . The operator j£ reads 

= ^EE^^M^*)-^-]- (5) 

R,a 5 

The current term of the operator proportional to A x (t) will 
not be used in this paper, and therefore it is pointless to give 
its form hereM^ 



B. Phonon modes 

In order to describe the effect of phonons in graphene we 
adopt the model developed by Woods and Mahan 47 and ex- 
tensively used by other authors i 48 i 49 ' 50 - 51 The potential energy 
of the model is a sum of two terms. The first is due to bond- 
stretching and reads 

y i = f E E t™^) - ^ <v a ] 2 , (6) 

R 6 

where u A (R) and ub(R + 8) represent the small displace- 
ments relatively to the equilibrium position of the carbon 
atoms in the sub-lattice A and sub-lattice B, respectively. If 
only the term (O is used a simple analytical expression for the 
eigen-modes is obtained. 52 The second term of the model is 
due to angle deformation and has the form 



V * = ^4 E {t^/ 2 + *0 ' ( UA ^ - UB ( R + + + «i) ■ (u A (R) - u B (R + d 2 ))} 2 

+ [(tfi/2 + S 3 ) ■ (u A (R) - u B (R + 8i)) + (S 3 /2 + i x ) • (u A (R) - u B (R + S 3 ))] 2 
+ [(5 2 /2 + S 3 ) ■ (u A (R) - u B (R + S 2 )) + (S 3 /2 + 8 2 ) ■ (u A (R) - u B (R + S 3 ))] 2 
+ [{St/2 + 8 2 ) ■ {u B {R + S 3 ) - u A {R - oi)) + (5 2 /2 + 8 1 ) ■ {u B {R + 8 3 ) - u A (R - a 2 ))f 
+ [{St/2 + 8 3 ) ■ {u B {R + 8 3 ) - u A {R - ai )) + (Si/2 + <M • {u B {R + 8 3 ) - u A {R))} 2 
+ [(8 2 /2 + 8 3 )-(u B (R + S 3 )~ UA (R-a 2 )) + (8 3 /2 + S 2 )-(u B (R + S 3 )-u A (R))] 2 } , (7) 

I 

where ax = <5i — ^3 and a 2 = 8 2 — 83, represent the vectors where Mq is the Carbon atom mass. The Lagrangian L = 
defining the unit cell. The kinetic energy has the form T — V\ — V2 leads to an eigenproblem of the form Mlo 2 w = 

T = Y,^a(R)} 2 + ^[ub(R)}\ (8) 
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Dw, where 



w = 




(9) 



and D is the dynamical matrix reading 



D 



/ X l F A C 

F* X 2 C B 

A* C* X x F 

V C* B* F X 2 



(10) 



with 



X 1 



Xo 



^-(3 + -/3 cos(V3Q y a) , 
o o 

45 3 

— /3 - -f3cos(V3Q y a) 

o o 

/3 cos(3Q x a/2) cos(V3Q y a/2) , 



(11) 



(12) 



A = -ae- l ^ a - ff + cos(V3Q J/ a/2)e ig » Q /^3) 

S = -^/3e-^ a - Qa + ^cos(V3g y a/2)e^fft) 

C = -iV3 f| - sin(V3Q y a/2)e iQ * a / 2 , (15) 

F = i/3^[sin(V3g y a)-2sin(V3g y /2)e 43Qa = a/2 ],(16) 
8 

(17) 

where we have redefined (3 as /3/a 2 and Q = (Q x , Q y ) is the 
momentum of the excitation. This model can be diagonalized 
numerically and its spectrum is represented in Fig. [TJ 




FIG. 1: (color online) Phonon spectrum of the honeycomb lattice, 
using a — 500 N/m and j3 = 10 N/m, such that the optical frequency 
at the T-point is or the order of 1600 cm~ 1 . Also shown are the spec- 
tra of the simplified model discussed in the text (dashed lines) and the 
effective model discussed in Appendix lAl (dashed-dotted lines). 

Although this model can easily be solved numerically, it is 
useful to derive from it a simple analytical model which helps 



in the analytical calculation of the electron-phonon problem. 
To do so, we follow Suzuura and Andc 48 and introduce two ef- 
fective models, one for the acoustic and the other for the opti- 
cal excitations, associated with the fields u = {ua + ub)/V% 
and v = (ua — ub)/V%, respectively. The effective Hamil- 
tonian for the acoustic modes has the form 



H = K ia 2 



Qx QxQy 
^xQy Qy 



-K 2 a J 



Ql + Q 
o 



o 

fx + Ql 



(18) 

with K x = 3a/8 and K 2 = 27a(3/{8a + 36/3/2). The eigen- 
modes of this effective Hamiltonian are 



M c lu 2 =K 2 (Qaf 



and 



M c u 2 = {K 2 + Kl){Qa) 2 , 
with polarization vectors 

£TA = i(—Qy, Qx)/Q, 

and 

CLA = i(Qx, Qy)/Q , 

respectively. The velocity of the modes is given by 



v TA = ay/K 2 /M c , 



and 



VLA 



= a^/{K x + K 2 )/M c 



(19) 



(20) 



(21) 



(22) 



(23) 



(24) 



The eigenmodes ( fT9l and (|20T > are represented in Fig. Uby 
dashed lines. They cannot be distinguished from the effective 
model described in Appendix lAl The effective Hamiltonian 
for the optical modes is given by 



H = K ( ; - ] - K 3 a i n n n2 



KaO, 



1 
1 



it Ql + Ql 



1 ( Qx QxQy 
\ QxQy Qy 



Ql + Q\ 



(25) 



with K = 3(a + 9/3/2), K 3 = (3a/ 8 - 27/3/8), and K 4 = 
27/3(a + 9/4/3)/ (8a + 36/3). The eigenmodes are 



M c to 2 =K -K 4 (Qa) 



and 



M c u? = K - (K 3 + K^iQa) 2 , 
with polarization vectors given by 

ETO = i(~Qy, Qx)/Q, 

and 

CLO = i(Qx, Qy)/Q , 

respectively. The modes d26b and d27b are represented in Fig. 
(T} by dashed lines. The dashed-dotted lines of the same 
Figure are derived from an effective model, discussed in Ap- 
pendixlAl 



(26) 
(27) 
(28) 
(29) 



4 



C. Electron-phonon interaction 



III. ELECTRONIC GREEN'S FUNCTION 



We now address the question of the electron-phonon inter- 
action. This comes about because the hopping t depends on 
the absolute distance between neighboring carbon atoms. We 
therefore have 



t[a + (u A (R) - u B (R + S)) ■ 8 /a] ~ t[a] 
1 dt[a] 
a da 



l -(u A (R)-u B {R + 8))-8. 



(30) 



Replacing (f30b in the Hamiltonian ([T]) and introducing the 
Fourier representation 



a tr (R) = -^-^e ih - R a a (k), 
V iv r , 



and 



(31) 



(32) 



with similar equations for b a (R) and ub(R), the electron- 
phonon interaction has form 



H, 



1 dt 1 



e—ph 



a da y/N c 



^2^2[ u a{Q) - u B (Q)e tQ ' s } ■ 8 



Q.k <y,S 

x [e ife ' 5 4(fe + Q)b a {k) + e- lk - s bl(k)a a (k - 0}) 

Since we are interested in the effect of the phonons with mo- 
mentum close to the r -point, the phase e lC *' s is expanded as 
e iQ S ^ \ — iQ ■ 8. Introducing the optical modes v, the 
electron-phonon interaction with the optical phonon modes 
has the form 



£ - ph ada^N c j-i ^ V McMQ) 

<- Q.k cr,i/,d V 



eAQ) ■ s 



x {Bl Qu + B Q ^[e lk - d al{k + Q)b a {k) 



— ifc-<5 > 



e 



>Uk)a a (k-Q)] 



(34) 



Introducing the acoustic modes u, the electron-phonon inter- 
action with the acoustic phonon modes has the form 



The electronic Green's function in the Dirac cone approxi- 
mation has the form 



<G 



iuj n — E(zw n ) —t<f)(k) 
—t<p*{k) iuj n — Tt(iui n ) 

[iu> n ~ E(£w n )][iw„ - T,(iu> n )] - t 2 \4>{k)\ 7 



(36) 



with <fi(k) — ^2ge lk S - The electronic self-energy shall be 
given by 

E(iw n ) = ?* m *{K, iu n ) + TP ho {K, iu n ) , (37) 

where E imp (K,iw„) = E unit (iw„) + Ec ost ;(-K", iuj n ) rep- 
resents the contribution due to mid-gap states (unitary scat- 
terers) as well as long-range Coulomb scatterers. The self- 
energy E pho (K,ioj n ) = Y, opt (K,iuj n ) + Y, ac (K,iuj n ) rep- 
resents the contributions due to optical and acoustic phonons. 
We note that the self-energies originating from the electron- 
phonon and Coulomb interaction are evaluated at the Dirac 
momentum K = 27r(l/3, \/3/9)/a. I n the following, we 
will discuss the Green's functions due to the the various con- 
tributions. 



A. Green's function with mid-gap states (unitary scatterers) 

The physical origin of mid-gap states in the spectrum of 
graphene is varied. Cracks, edges, vacancies 53 are all possible 
sources for mid-gap states. From an analytical point of view, 
these types of impurities (scatterers) are easily modeled by 
considering the effect of vacancies. We stress, however, that 
this route is chosen due to its analytic simplicity. 

The effect of mid-gap states on the conductivity of graphene 
was first considered by Peres et al. 20 for the case a half-filled 
system. Considering the effect of a local scattering potential 
of intensity eo, the Green's function has the form of Eq. 
with the retarded self-energy 



n;e 1 - eoF{uj) - ie 7ri?(w) 



h [1 - e F M] 2 + [e 7ri?H] 2 
where the functions F(uj) and R(u>) are defined by 

£ G(k, u + i0+) = F(lo) - ittR(lo) . 



(38) 



(39) 



rrac 
e—ph 



dt 1 



EE 



a da VNc jri V 4Mcrw„(Q) 

v c Q.k a,v,o V 

x (Q • 8){Bl Q v + B Q .„){e ik S aUk + Q)b a (k) 
+ e- lk - s b\(k)a a (k-Q)\. (35) 

In both Eqs. d34l i and ( f35T > the Bq v (Bq v ) operators are de- 
struction (creation) operators of phonons of momentum Q and 
polarization v. 



Mid-gap states are obtained by making the limit eo — > oo, 
which resembles the unitary limit. Clearly, R(u>) is the den- 
(e !/ (Q) ■ 8) sity of states per spin per unit cell. Let us write S^* it (ui) as a 



sum of real and imaginary parts, E£^ t (w) = S'(w) + £E"(cj) 
(note that E" > 0). The functions F(u) and R(to) are deter- 
mined self-consistently through the numerical solution of the 
following set of equations: 



F(w) - mR(w) 



1 



:[7iE"(T-i#/2) 



-(fej-nS')(*/2-iT)], (40) 
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with T and W given by 

hH' — coH 



arctan ■ 



a=±l 



arctan ■ 



KS' -Luh + aD 



and 



(41) 



(42) 



In Figure HI we compare the density of states computed 
using the coherent potential approximation (CPA) equations 
with that obtained from a numerical exact method. 54 It is clear 
that the CPA captures the formation of mid-gap states in a 
quantitative way. The main difference is the presence of a 
peak at zero energy in the exact density of states, whose mea- 
sure is quantitatively negligible. 
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FIG. 2: (color online) Density of states of graphene in the presence 
of mid-gap states. The CPA calculation is compared with a numerical 
exact method. The concentration of impurities is n% = 0.005 and 
m = 0.01. Here and in the following figures we use t — 3 eV and a 
cutoff energy of D — 7 eV. 

In Figure [3] we depict the self-energy calculation using the 
CPA equations, for two different values of the impurity con- 
centration. It is clear that Tj"(u>) increases close to zero en- 
ergy leading to a broadening of the electronic states close to 
the Dirac point. 



B. Green's function with Coulomb impurities 

It has been argued that charged impurities are crucial to 
understand the transport properties of graphene on top of a 
silicon oxide substrate i 31,55 ' 56 In what follows, we compute 
the electronic self-energy due to charge impurities, using sec- 
ond order perturbation theory in the scattering potential. Elec- 
tronic scattering from an impurity of charge Ze leads to a term 
in the Hamiltonian of the form 



V = - 



Ze 1 



R 2 



[ai(R)a a (R) + bi(R)b a {R)]. (43) 



n =0.001 

n =0.0005 0.04 




-6 -4-2 1 



FIG. 3: (color online) CPA calculation of the self-energy hT^ it (u>) 
for two impurity concentrations, m = 0.001 and m — 0.0005. The 
left panel shows the imaginary part and the right one shows the real 
part of the self-energy. 



In momentum space V reads 

V=^~ J2 Vo(QM(pW(p + q) + bUp)b a (p + q)], 

(44) 



where Vq (q) reads 



Vo(q) = -J2 



R 



Ze 2 eJ R i 
Vd 2 + R 2 



(45) 



With <G°(k, iu! n ) the bare and G(k, p, iui n ) the full Green's 
functions, the Dyson equation due to one Coulomb impurity 
reads 

G(k,p,w n ) = S k .p + G°(k,iuj n ) 

c k' 

If we consider a finite density per unit cell, nf, and inco- 
herent scattering between impurities, the second-order self- 
energy is given by 



h?N, 



J2v 2 (k-p)<& (p,iu; n ), (47) 



where a term of the form nfV(0)/H was absorbed in the 
chemical potential, since it corresponds to an energy shift 
only. Note that we have replaced Vo(q) by V(q), which cor- 
responds to include the effect of electronic screening in the 
calculation. The form of V(q) is (in S.I. units) 33 



V(q) 



Ze 2 e~ qd 
2e eA c q + 7 



(48) 



where e = 3.9 is the Silicon Oxide relative permittivity, d is 
the distance from the charge to the graphene plane, and 7 is 
given by 



7 



(49) 
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where p(p) is the self-consistent density of states as computed 
from the CPA calculation ( A c = 3\/3a 2 /2 is the area of the 
unit cell). 

The self-energy d47| i is dependent both on the momentum 
k and on the frequency. However, we are interested on the ef- 
fect of the self-energy for momentum close to the Dirac point. 
Within this approximation the imaginary part of the retarded 
self-energy becomes diagonal and momentum independent, 
reading (d ~ 0) 



Coul 



Z 2 e i n? , , /2|M 

\hu\ I -4r— - +7 



4A2 e 2 e y^a 1 



3ta 



(50) 



The self-energy contribution due to Coulomb impurities is 
the most relevant one in order to fit the experimental data of 
Ref. B26I1 and is shown on the right hand side of Fig. [6] 



C. Green's function with phonons 

Following the same procedure as in the previous subsec- 
tion, the self-energy due to optical phonons within first order 
perturbation theory is given by 



2 \da J hM c to N c 



Q 



l - J2 D°(Q, iu m )G°(K - Q, icj n - w m ) . (51) 



The analytical form of the self-energy due to acoustic 
phonons reads 

iuj) 9 ( dt Y 1 V {Qa? 
1 ' n) ~ 16 \da) N c 4HMcUu(Q) 

x JfiJ2 iv m)G°\K - Q, uu n - iv m ) . (52) 



The unperturbed Green's functions have the form 

2w„(Q) 



D"{Q,iv m ) 



(iv m y - K(Q)] 2 ' 



(53) 



approximates by uo v (Q) ~ ujq, the calculations are simple to 
do and the result for the imaginary part of the self-energy due 
to optical phonons is 



A a A a - s/3D 
2 A n arctan — 2A n arctan 

[ns(7la;o) + (1 — a)/2 + anp^hw + aHujQ — //)] , (55) 



with s = sign A, 



A Q = So? + a/iwo — fiS (a; 
7" = hY;"(uj + auj ) , 



CtLOo) 



(56) 
(57) 



and tib(x) and uf(x) the Bose and Fermi functions, respec- 
tively. 

The self-energy due to optical phonons, computed using 
the disordered electronic Green's function, is compared with 
the same quantity computed using the bare electronic Green's 
function in the central panel of Fig. [4] It is clear that the 
imaginary part of the self-energy has a larger value when the 
disordered Green's function is used. However, in the region 
of frequencies // — Hluq < Hlj < Hujq + //, at T = 0, the 
imaginary part coming from the optical phonons is zero, both 
when one uses the bare and the disordered Green's functions. 
This is due to the arguments of the Bose and Fermi functions. 



11=0.29 eV, T=45K 



and 



m (eVj 




G°(k,iLO n ) 



{iLu n y-t^{k)\ 2 ' 



(54) 



and uiq — y/Ko/Mc- The Matsubara summation over the 
frequency v m is done using standard methods. 



1. Effect of Disorder: Fermionic Propagator 

If we include the effect of disorder, the unperturbed Green's 
function G°(K — Q, ioj n — iv m ) should be replaced by the 
dressed Green's function due to the impurities. For the (rele- 
vant) case of optical phonons, where the phonon dispersion is 



FIG. 4: (color online) Electronic self-energy due to unitary scatterers 
and acoustic and optical phonons. Only the imaginary part is repre- 
sented. The impurity concentration is m = 0.0002. The self-energy 
due to optical phonons depends on temperature and on the chemi- 
cal potential p. We have chosen p = 0.29 eV and T = 45 K. In 
the central panel we show the calculation for self-energy due to op- 
tical phonons both using the bare electronic Green's function (solid 
line) and the Green's function with mid-gap states (dashed line). The 
acoustic self-energy is independent of the chemical potential and is 
represented for a temperature of T — 45 K. 

In Fig. |4] we depicted the self energy of the short ranged 
impurities together with those due to acoustic and optical 
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phonons. It clear that the effect of acoustic phonons is negli- 
gible at low energies. The self-energy due to optical phonons 
depends on the chemical potential, and is represented for a 
gate voltage of V g = 71 V O = 0.29 eV). 



2. Effect of Disorder: Phononic Propagator 

To be consistent, also the phonon propagator has to 
be dressed due to its interaction with impurities. The 
phonon propagator shall be renormalized within the RPA- 
approximation, i.e., 



D 



RPA 



(iuj n ) 



2w„ 



(iu n ) 2 - (w„) 2 - 2w„n i ,(jw n ) 



(58) 



where the first order of the phononic self-energy H„(u> n ) is 
proportional to the polarization defined as 



i r h/3 

P«K) = lim — / dTe^ T (T T p(q,T)p(-q,0)) 

q^O A s J Q 



(59) 



with p(—q) denoting the density operator. Explicitly, we get 
for the imaginary part of the retarded phononic self-energy 



dt\ 18 



Muj v \daj y/3irt 2 J 2vr 
x [np(u!ih — p) ~ np(u!ih + uh — p)] , (60) 



where the dimensional function Q(u>i,u>) is given in Ap- 
pendixlBl 

In Fig. |5] the self-energy due to electron-phonon scatter- 
ing (left) and the conductivity (right) are shown as they result 
using the bare (dashed) and dressed (full) phonon propagator. 
Since the effect is hardly appreciable, we show also the re- 
sults where the phononic self-energy has been multiplied by a 
factor 10 (dotted-dashed). The renormalization of the phonon 
propagator due to disorder is thus negligible and the results of 
the following section will be obtained using the bare phonon 
propagator. 



IV. THE DC AND AC CONDUCTIVITY 

In this section, we discuss the transport properties of 
graphene due to the various sources of one-particle scattering. 
This is done within the Kubo formalism. 



00013, N=0.00004, T=45K, V=71 Volts 




— bare phonon propagator 

— dressed by n((o) 

— dressed by JfMI(co) 




2000 4000 6000 8000 

-I 

(o in cm 



FIG. 5: (color online) Electronic self-energy due to electron-optical 
phonon interaction 3Sop* (left) and optical conductivity a (right) for 
various phonon propagators. The electronic propagator is dressed 
by the self-energies due to Coulomb (nc = 0.00013) and impurity 
scattering (m = 0.00004). The temperature is T — 45K and the 
applied gate voltage V = 71Volts. 



with A s = N C A C the area of the sample, and A c — 3\/3a 2 /2 
the area of the unit cell, from which it follows that 



$l(Txx{u>) = D8{lo) + 



hujA s 



and 



$sa xx (u>) = - 



A s uj fuvA s 
where D is the charge stiffness which reads 

MA xx (lj + i0+) 



D 



-IT : 7T- 



A. 



hA s 



(62) 



(63) 



(64) 



The function A^a; + i0 + ) is obtained from the Matsubara 
current-current correlation function, defined as 



A xx (iu! n ) = I dre 1 ' 
la 



<Trjf(T)jf(0)> 



(65) 



The calculation of the conductivity amounts to the determina- 
tion of the current-current correlation function. 



B. The real part of the DC conductivity 



A. The Kubo formula 

The Kubo formula for the conductivity is given by 

A xx (Lu + iO+) 



< if > 



iA s (uj + i0+) ihA s (u + i0+) ' 



(61) 



The real part of the DC conductivity is given by 

=~l deK{e) df{€ -^ , (66) 
nn J oe 

where f(x) is the Fermi function and K(e) is a dimensionless 
function that depends on the full self-energy. In the limit of 
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zero temperature the derivative of the Fermi function tends to with K(fi) given by 
a delta-function and the conductivity is given by 

5R<t(m) = ^j-K(fi) , (67) 



_,, , 1 / D(sA + D)I a A A A + sD 

K ^ = W a 2^ { (A + sDY + Zg + ^ axctan - - A arctan 

7 / A A + D 

+ —2- ( 2 arctan arctan ^4 — D/ a — arctan 

2AV I a I a 



and 

4 = ?£y£(/i) ,A = (i- mtE{n) (68) 

and D the cutoff bandwidth. 

In Fig. [6] we plot ctdc(aO as function of the gate voltage, 
Vg, considering the effects of both charged impurities, mid- 
gap states and acoustic phonons (the optical phonons do not 
contribute to cfdc{^))- It is clear that the linear behavior of 
cue (a*) as function of V g is recovered. The theoretical curves 
are plotted together with the experimental data of Ref. B26I1 . 
After having fitted the experimental curves, it is clear that the 
most dominant contribution to the DC conductivity is coming 
from Coulomb impurities. 




V g (volt) to (eV) 

FIG. 6: (color online) Left panel: Experimental data (circles) from 
Ref. 12611 . conductivity with the full self-energy (solid line) due 
to phonons, mid-gap states, and charged impurities, conductivity 
(dashed-line) with self-energy due to mid-gap states and charged im- 
purities, only. The parameters are T — 45 K, m — 4 x 10~ 5 , 
nf = 7.5 x 10 -5 , and d = 0. The dashed-dotted line is the same 
as the solid line but with nf = 1.3 x 10~ 4 . Right panel: Imaginary 
part of the self-energy, 7iSrE£f ^ (K, to), due to charged impurities, 
the most dominant contribution. 

Let us now further discuss the conductivity as function of 
the gate voltage V g , which relates to the chemical potential 
as V g oc fj 2 . In Figure [7] we show a(/j,) as function of V g 
considering both charged impurities and short range scatter- 
ers (left panel), having the same impurity concentration. We 



note that the conductivity, albeit mostly controlled by charged 
impurities still has finger prints of the finite eo scatterers, due 
to the asymmetry between the hole (negative V g ) and particle 
(positive Vg) branches. The conductivity follows closely the 
relation er(/i) oc V g , except close to the Dirac point where is 
value its controlled by the short range scattering. 

If we suppress the scattering due to charged impurities, 
which should be the case in suspended graphene, only the 
scattering due to short range scatterers survive. In this case the 
right panel of Fig. [7] shows that there is a strong asymmetry 
between the hole and the particle branches of the conductivity 
curve even for a value of eo as large as 100 eV. Moreover, the 
smaller the value of eq the larger is the asymmetry in the con- 
ductivity curve. We note that asymmetric conductivity curves 
were recently observed. 57 On the contrary, if the experimental 
data shows particle-hole symmetry of <r(/z) around the Dirac 
point, then the dominant source of scattering is coming from 
very strong short range potentials, i.e., scatterers that are in 
the unitary limit. 



nf=n F 0.0001 nf=0, n F 0.0001 




V g (volt) V g (volt) 



FIG. 7: (color online) Left: Conductivity <t(/i), in units of ao = 
Tve 2 /(2h), considering both short range and charged scatterers. 
Righ: Conductivity cr(u) when the influence of charge scatterers is 
removed. 
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C. The real part of the AC conductivity 

The finite frequency part of the conductivity is given by 

3?cr(w) = ^— [ — 8(wi,w)[/(fo 1 -/i)-/(foi+/KJ-/i)] 
irn J uj 

(69) 

where f{x) is the Fermi function and 0(ci>i,w) is a dimen- 
sionless function given in Appendix iBl As can be see from 
this Appendix, the function 9(a;i, uj) depends on the self en- 
ergy, which is due both to impurities and phonons. 



1. Optical conductivity without Coulomb scatterers 

We will first discuss the optical conductivity without 
Coulomb scatterers since they should not be present in sus- 
pended graphene. In Figure [8] we plot the conductivity of a 



graphene plane in units of ctq 



7T JX 



I h. The calculation is 



made at two different temperatures, T — 45 K and T = 300 
K, and for a density of impurities m = 0.0004. The calcu- 
lation compares the conductivity with and without the effect 
of the phonons. The main effect induced by short-ranged im- 
purities is the existence of a finite light-absorption in the fre- 
quency range < hui < 2/i. The optical phonons increase 
the absorption in this frequency range. We have checked that 
the effect of the acoustic phonons is negligible. The optical 
phonons also induce a conductivity larger than oo for frequen- 
cies above u> = 2/i. The effect is more pronounced at low 
temperatures. 



f? = 0.0004, 7=45 K 
with phonons 




5000 10000 
co (Mem) 




side). There is more absorption for frequencies in the re- 
gion < hw < 2/i the larger the impurity concentration is. 
This is because the number of mid-gap states is proportional 
For larger gate-voltage a plateau is reached at fre- 



to 



20 



quencies in < hw < 2/i. For small gate voltages, V g , the 
absorption in the region < hw < 2/i is larger than for larger 
V a . 



V=71 volt, T=45K 



n =0.0004, T=45K 




5000 10000 
CO (1/cm) 



5000 10000 
co {1/cm) 



FIG. 9: (color online) Optical conductivity of graphene at temper- 
ature T =45 K. The left hand side shows curves for various short- 
ranged impurity densities rii . The right hand side shows curves for 
various gate voltage V g . The vertical, dashed lines mark twice the 
chemical potential. 



0.00013, n = 0.00004, T = 45 K 




5000 10000 
CO (1/cm) 



2000 4000 6000 

co (1/cm) 



FIG. 8: (color online) Optical conductivity of graphene at two dif- 
ferent temperatures T —45, 300 K. Each panel has two curves, the 
solid curve is the conductivity with only short-ranged impurities, the 
dashed curve is the conductivity with both short-ranged impurities 
and phonons, with the self-energy due to phonons computed with the 
full Green's function. The chemical potential is that associated with 
a gate voltage of V g — 71 V. The vertical, dashed line marks twice 
this value. 

In Figure [9] we again plot the conductivity of a graphene 
plane in units of <tq — ^e 2 /h and at temperature T = 45 
K. This time we compare different impurity densities rii (left 
hand side) and gate voltages/chemical potentials (right hand 



FIG. 10: (color online) Real part of the optical conductivity including 
the effect of phonons, midgap states and charged impurities. The 
parameters are T = 45 K, m = 4.0 x 10~ 5 , and nf = 1.3 x 10~ 4 . 
The dashed vertical lines correspond to hui = 2/i and the shorter 
dotted-dashed ones to hui = 2hu>o + 2/i. 



2. Optical conductivity with Coulomb scatterers 

We will now discuss the optical conductivity with Coulomb 
scatterers which are generally present in graphene on a 
substrate^. In Figure [10] we plot the conductivity of a 
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nf = 0.00013, n, = 0.00004, T=45K 




2000 4000 6000 8000 



oo (1/cm) 

FIG. 1 1 : (color online) Real part of the optical conductivity including 
the effect of phonons, mid-gap states and charged impurities, com- 
pared with the experimental data from Ref. I26TI . The parameters are 
T = 45 K, n, = 4.0 x 10~ 5 , nf = 1.3 x 10" 4 , and d = 0. 

graphene plane in units of uq = f e 2 //i, including the ef- 
fect of midgap states, charged impurities and phonons. Again, 
the main feature is that the conductivity is finite in the range 
< huj < 2/i and increases as the gate voltage decreases. 
We choose the concentration of midgap states in Fig. [TOlto 
be one order of magnitude smaller than the one of Coulomb 
scatterers^, and therefore the conductivity is mainly con- 
trolled by phonons and charged impurities. 

Another feature of the curves in Fig. [I0]is the large broad- 
ening of the interband transition edge at tiui = 2// (indicated 
by vertical dashed lines). Note that this broadening is not due 
to temperature but to charged impurities, instead. In fact, the 
broadening for all values of V g is larger when the conductiv- 
ity is controlled by charged impurities rather than by midgap 
states. 

The coupling to phonons produces a feature centered at 
2hidQ + 2/i where hujQ denotes the LO-phonon energy corre- 
sponding to the wave number 1600cm -1 (indicated by shorter 
dotted-dashed vertical lines). For gate voltages with /i < Hu>o 
there appears a similar feature at ITiloq — 2/i which is not 
washed out by disorder. The optical phonons thus induce a 
conductivity larger than <7n around these frequencies. This ef- 
fect is washed out at larger temperatures and a < ctq. 

All these effects are consistent with the recent infrared mea- 
surements of graphene on a S1O2 substrate 26 which are shown 
as dashed lines in Fig. Q~T] Notice that there is only one fit- 
ting parameter involved which is adjusted by the conductivity 
curve at zero chemical potential. Whereas for low gate voltage 
the agreement is good, there is considerable weight missing 
for higher gate voltage (V > 20 Volts). Nevertheless, all the- 
oretical lines predict lower conductivity than the experimen- 
tal measurements, such that our model with only one fitting 
parameter is consistent. Since we have included in our cal- 
culation all possible one-particle scattering mechanisms, the 
missing weight at large gate voltages could be attributed to 
electron-electron interactions in graphene, which become im- 
portant at this electronic densities. 



V. DISCUSSION AND CONCLUSIONS 

In this paper we have computed the optical conductivity of 
graphene at finite chemical potential, generalizing the results 
of Ref. 112011 . The calculation includes both the effect of dis- 
order (mid-gap states and charged impurities, which have a 
different signature in the DC conductivity 58 ) and the effect of 
phonons (optical and acoustic). It is shown that at low tem- 
peratures the effect of acoustic phonons in negligible, since it 
induces an imaginary part of the electrons' self-energy that is 
much smaller than the imaginary part induced by the impuri- 
ties. For a discussion based on the Boltzmann equation, see 
Ref.lH 

The imaginary part induced by optical phonons is of the or- 
der of the imaginary part induced by impurities. Still, optical 
phonons are only important in the calculation of the real part 
of a{ui), they play no role in the calculation of ctdcC/- 4 ) m tne 
temperature range T G [0, 300] K. The self-energies due to 
mid-gap states and due to optical phonons are in a sense com- 
plementary, since the imaginary part coming from impurities 
is large at the Dirac point, whereas the imaginary part com- 
ing from the optical phonons increase linearly with the energy 
away from the Dirac point. The imaginary part of the self en- 
ergy due to charged impurities is a non-monotonous function 
of the energy, growing first linearly but changing to a decaying 
behavior of the form l/ui, for large energies. 

In Section [IV] we have discussed the different scatter- 
ing mechanisms separately, since for suspended graphene 
Coulomb scatterers will be absent. Since vacancies, corre- 
sponding to an infinite potential where particle-hole symmetry 
is restored, are unlikely, we model the short-ranged potentials 
due to cracks, ripples, etc. by large, but finite short ranged po- 
tential. For the DC conductivity, this leads to an asymmetry 
of the hole- and electron-doped regime. 

The most notable effect on the optical conductivity of both 
the optical phonons and the impurities is the induction of a 
finite energy absorption in the energy range < hut < 2/i, 
a region where the clean theory predicts a negligible absorp- 
tion. The optical phonons also induce a conductivity larger 
than (To around fko ~ 2/i. It is interesting to note that for 
frequencies away from the Dirac point the imaginary part of 
the self-energy due to optical phonons is linear in frequency, 
a behavior similar to that due to electron-electron interactions 
in graphene. 

It is clear from Fig. QT|that m general the calculated ab- 
sorption in the range < hio < 2/i is not as large as the 
experimental one. It is also noticeable that for small gate volt- 
ages there is a reasonable fit of both the absorption and of 
the broadening of the of the step around 2/i. For large val- 
ues of the gate voltage the calculated absorption is smaller 
than the measured one. This suggests that some additional 
scattering mechanism is missing in the calculation of the op- 
tical conductivity. The missing mechanism has to be more 
effective at large gate voltage. A possibility are plasmons of 
the type found in Ref. 15911 . Another possibility is that water 
molecules which are especially active in the infrared regime, 
are contributing to the missing weight. 

In synthesis we have provided a complete and self- 
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consistent description of the optical conductivity of graphene 
on a substrate (including Coulomb scatterers) and suspended 
(without Coulomb scatterers). Our results are in qualitative 
agreement with the experimental results. To meet a quantita- 
tive agreement further research (also on suspended graphene) 
is necessary, but water molecules underneath the graphene 
sheet are likely to account for the missing weight. 



eigenvalue problem can be put in the form 



B\U + B 2 v = MctjJ u 



Bin + B 3 v = M c uj 2 v 



(A5) 
(A6) 



Let us first look at the acoustic modes. They corresponds to 
oj ~ 0. In this case we can write 
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-u , 



3(a + 9/3/2) 
from which an eigenvalue equation for u follows: 



BoBn n 

B\U ; 2 ; X M = MU?U . 

3(a + 9/3/2) 



(A7) 



(A8) 



In the case of the optical modes one has Alio 2 ~ 3a + 27/3/2, 
from which we can write 



APPENDIX A: AN APPROXIMATE PHONON MODEL 



Bo 



u 



3(a + 9/3/2) 



(A9) 



We start by expanding the matrix elements of the dynamical 
matrix ( [Tol l up to second order in the momentum Q. We then 
introduce acoustic, u = (ua + ub)/V2, an d optical, v = 
(ua — ub) I v2, modes. This procedure leads to the following 
eigenvalue problem 



M c uj 



Bi B 2 
Bl Bi 



u 



(Al) 



which leads to 



B\B 2 



3(a + 9/3/2) 



v = Muj 2 v . (A10) 



The dashed-dotted lines in Fig. Q]are the eigenvalues of Eqs. 
and (lATOl . 



with 



' 32 I kxQxQy y 2 



AaQ x Q y 



B 2 = 



(A2) 



with yi = 3(2a + 3/3)Q 2 + (2a + 9f3)Q 2 y and y 2 = (2a 
9/3) Q 2 X + 3 (2a + 3/3) Q 2 Further, we have 



(A3) 



and 



o _ 3^ ( Vs (36/3 - Aa)Q x Q y 

3 32 \ {mp-Aa)Q x Q v V i 

' (A4) 

with y 3 = 32a + 144/3 - (6a + 9f3)Q 2 x - (2a + 45/3)Q 2 ( and 
2/3 = 32a + 144/3 - (6a + 9/3)Q 2 _ (2a + Ah(3)Q 2 x . The 



APPENDIX B: THE FUNCTION 6 (wi, u) 

In Eq. (|69b the dimensionless 0(u;i, oj) function was intro- 
duced. Let us define 



A = fuvi — KY/(u>i) , 
B = fvjji + huo — fiS (u>i + oj) , 
I a = KE"(wi) , 
h = hT,"(uj 1 +uj), 
Di = 2[{A-B) 2 + {I a -I b f] x 

[(A-B) 2 + (I a + I b ) 2 ], 
D 2 = 2[(A + B) 2 + (I a -I b ) 2 } x 
[(A + Bf + (4 + I b ) 2 ] . 



(Bl) 
(B2) 
(B3) 
(B4) 

(B5) 

(B6) 



The function Q(uj%, oj) is expressed in terms of the above aux- 
iliary functions as follows 
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9(>i, uj) = -^-h(A 3 -2A 2 B -2BI 2 + A(B 2 + I 2 + I 2 )) ^ (arctan — - arctan A + aD \ 

D l a = ±1 \ la la J 

2 (A 
+ jyh{B 3 - 2B 2 A - 2All + B{A 2 + I 2 + I' 2 )) ^ (arctan — - arctan 

O . , / A 

+ ft 



"-^(A 3 + 2A 2 B + 2B/ 2 + ^(B 2 + / 2 + J 2 )) V arctan — - arctan - 

J 2 _±Y n V i 





la 


.4 


+ aD 




h 


.4 


+ aD 




la 


.4 


+ aD 



a=±l 

2 x / A 

—I a (B 3 + 2B 2 A + 2AIj ! + B{A 2 +ll+I 2 )) J2 (arctan — - arctan , 

a=±l 

±,±\tt.,a2 „2, T 2 T 2,^(^A 2 + Il ,„„ (A — aD) 2 + I 2 

ft ft 



w - B ' + < - J > > h - 108 <b 1 7u? + n ) ■ <B7) 
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